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Abstract: We discuss collider signatures of the "natural supersymmetry" scenario with 
baryon-number violating R-parity violation. We argue that this is one of the few remaining 
viable incarnations of weak scale supersymmetry consistent with full electroweak natural- 
ness. We show that this intriguing and challenging scenario contains distinctive LHC 
signals, resonances of hard jets in conjunction with relatively soft leptons and missing en- 
ergy, which are easily overlooked by existing LHC searches. We propose novel strategies 
for distinguishing these signals above background, and estimate their potential reach at 
the 8 TeV LHC. We show that other multi-lepton signals of this scenario can be seen 
by currently existing searches with increased statistics, but these opportunities are more 
spectrum-dependent. 
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1 Introduction 

The scenario of weak scale supersymmetry (SUSY) with light third generation superpart- 
ners ("natural SUSY" or "effective SUSY") was put forward many years ago [1, 2] as a way 
of ameliorating the SUSY flavor and CP problems, while maintaining electroweak (EW) 
naturalness. The central observation from the hierarchy problem viewpoint is that, in the 
Standard Model (SM), the radiative corrections which destabilize the Higgs potential are 
dominated by the heavy top quark. EW stability requires SUSY cancellations from stops 
of a few hundred GeV. However, from a purely bottom-up viewpoint it is attractive to 
consider the other superpartners to be substantially heavier, since they are primarily im- 
plicated in the SUSY flavor and CP problems. The theoretical challenge is to then identify 
high-energy models and mechanisms which, at least approximately, lead to this kind of 
split superspectrum (see [3-9] for a partial list). 

There has been a resurgence of interest in this scenario in the last year, mainly moti- 
vated by the null results thus far of LHC searches for physics beyond the SM. In light of 
the current LHC bounds on R-parity conserving SUSY, it is difficult to envision any other 
viable version of SUSY which is consistent with full electroweak naturalness (that is, ab- 
sence of EW fine-tuning). By contrast, it has been recently shown that natural SUSY easily 
evades the most stringent LHC constraints with integrated luminosity C ~ 1 fb _1 [10-13]. 
Later dedicated searches for R-parity conserving natural supersymmetry have appeared, 
better constraining particular spectra with [14] or without [15] light gluinos, however the 
natural parameter space is still quite open. 

Natural SUSY has an even wider significance, in that it beautifully illustrates the 
general theme of how a "top-partner" can algebraically cancel destabilizing top-quark ra- 
diative corrections to the Higgs potential. This relates it to the theory and phenomenology 
of fermionic top-partners [16], appearing in non-supersymmetric Little Higgs (see for re- 
view [17, 18] )and Twin-Higgs models [19, 20]. In this sense, light stop searches fit into the 
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broader program of testing whether any top-partner is helping to stabilize the weak scale. 
The LHC is the first experiment in history that can test naturalness on such a broad front, 
in a relatively comprehensive and well-defined way Either a discovery of such top-partners 
at the natural scale of a few hundred Ge V, or even their exclusion to high confidence, would 
constitute a significant scientific finding. In previous work [11], we have argued that for this 
grand and challenging experimental undertaking, one should free the natural SUSY setting 
from too much UV prejudice and anticipation, lest this lead to overlooking experimental 
opportunities now and because UV considerations have not led to any sharp no-go "theo- 
rem". We have carved out a simple theoretical framework that facilitates this. The present 
paper will discuss important but widely overlooked signals that follow straightforwardly 
from this perspective, as well as the combination of methods that can separate them from 
SM background. 

In particular, although most experimental searches for natural supersymmetry have 
concentrated on the R-parity conserving case, we consider here the case of (baryon-number 
violating) R-parity violation (RPV). The plausibility and attractiveness of this scenario was 
argued in Ref. [11] from a number of viewpoints and considerations (also see [21] and refer- 
ences therein for recent models of spontaneous RPV) . The spectrum of RPV and R-parity 
conserving natural SUSY can be quite similar. If the theory is completely natural we expect 
both species of stops and at least the left-handed sbottom with masses of order ~ 400 GeV 
or lighter. Gluinos can be naturally twice as heavy if they are Majorana fermions, and 
even heavier if they are (part of) Dirac states [22]. That is, we cannot guarantee the gluino 
to be experimentally accessible in the near future. However, if we are lucky and Majorana 
gluinos are light enough in RPV natural SUSY, they can produce spectacular signatures 
in same-sign dileptons [23]. Here, we assume more minimally that the stops and sbottom 
mandated by naturalness are the only accessible colored superpartners. 

RPV is distinct from R-symmetry conservation because phenomenological viability 
does not require a neutral superpartner to be at the bottom of the SUSY spectrum, since 
superpartners are allowed to decay into SM particles. Therefore EW gauginos can easily be 
heavier than the stops and the sbottom, having no significant impact on the phenomenology. 

Non-minimal Higgs degrees of freedom are subtler. Higgsinos are usually assumed 
to acquire mass from the same "/ii-term" that also contributes to Higgs scalar potential. 
Higgs naturalness then requires Higgsinos not much heavier than ~ 200 GeV. (However, 
see Ref. [11] for a bottom-up description in which Higgsinos can be much heavier.) We 
will show in Sec. 2 that light Higgsinos can remain relatively well-hidden in the RPV 
context. On the other hand, extra Higgs boson degrees of freedom of SUSY can be heavier 
without compromising naturalness. We therefore only keep the SM Higgs scalar in our 
study of stop/sbottom phenomenology. Finally, there are by now stringent bounds on 
the SM Higgs mass, and even tentative hints of its presence at ~ 125 GeV. Theoretically 
accommodating the Higgs mass in high- energy models has been an increasing challenge 
ever since LEP2, but there are certainly interesting ideas for doing this. By contrast, a 125 
GeV Higgs mass is straightforwardly accommodated within our bottom-up natural SUSY 
framework, deferring the full UV description of physics lying outside experimental reach. 
In particular, we do not restrict stop/sbottom masses by their radiative contributions to 
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the physical Higgs mass, since there may well be other contributions from unknown heavier 
sources. 

Collider signatures of RPV SUSY are largely dictated by the detailed structure of RPV 
interactions, which cannot be anarchical (for a review see [24]). Either baryon number or 
lepton number should be conserved to avoid prompt single proton decay. While lepton 
number violation (LNV) is interesting by itself and deserves more study, it has already 
meaningful constraints from the LHC, since it mostly leads to leptons and taus in the final 
states, which are relatively easy to spot. Baryon number violation (BNV) is experimentally 
more challenging than LNV, resulting in jetty final states and suffering from enormous QCD 
and tt backgrounds. We will focus on this challenging scenario of two stops and a sbottom 
at the bottom of a BNV SUSY spectrum. Not only does this fill a gap in SUSY searches, 
but it also shares several features with, and insights into, other top-partner searches. This 
spectrum was considered earlier in Ref. [25] in the context of the CDF ll Wjj anomaly" [26] . 
Although the anomaly was later refuted by DO [27], as well as by a similar CMS search [28], 
this paper was an important step in understanding of collider signatures of the minimal 
spectrum. Here, we will elaborate on several points briefly touched on in Ref [25], broaden 
the motivations and scope, and detail new search strategies. 

A particular difference with Ref [25] is that we will assume that BNV is governed by 
small couplings and therefore we will neglect single-resonance production of superpartners, 
concentrating on pair-production. Such smaller couplings make the theory more straight- 
forwardly safe from low-energy precision data. We will argue that it is theoretically very 
plausible that a stop is the lightest superpartner, which can decay into a pair of jets. 
Needless to say, by itself this is an extremely challenging signature, given relatively small 
production cross sections and absence of any "interesting" features in the event, e.g. lep- 
tons or missing transverse energy (MET). However, one can take advantage of production 
of the heavier sbottom and stop which further cascade decay into the lightest stop, emitting 
W, Z and/or higgs (on- or off-shell) along the way. These events are more promising, be- 
cause they can potentially contain leptons and MET. Nonetheless, existing cut-and-count 
searches are not optimized for signatures like this and generally overlook them. They do 
not take advantage of the most important qualities of these events: hard jets reconstruct- 
ing a pair of resonances, in conjunction with leptons and/or MET which is relatively soft 
compared to the top quark background. 

We will substantiate these claims, and use them to craft a search strategy for the 
most promising and robust of these cascades. We will show that the backgrounds are 
under control and y/s = 8 TeV LHC can have a good reach for these events. We will also 
discuss other channels, which can be promising, but where the backgrounds are not easy to 
estimate with theoretical tools. The alternative possibilities for the lightest superpartner, 
a sbottom or Higgsino, are also plausible but even more phenomenologically challenging, 
and we defer their consideration from this paper. 

The paper is organized as follows. In the next section we review the BNV RPV 
natural SUSY scenario and reduce it to its most LHC-relevant features, thereby arriving 
at a useful "simplified model". In section 3 we study the various "charged current" (W) 
channels and relevant backgrounds, and roughly estimate which of these channels is viable. 
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In section 4 we perform explicit simulations of signal and background in the most promising 
of these channels and discuss cuts (which are quite different from standard SUSY searches) 
in greater detail. In section 5 we briefly discuss "neutral current" (Z, h) cascade decays 
between stops. While getting a substantial number of these events is only possible in a 
subset of stop/sbottom/higgsino spectra, they can be quite spectacular, and indeed the 
multi-lepton CMS search [29] already has an appreciable sensitivity. They remain an 
exciting discovery channel for the future with more statistics. Finally in section 6 we 
conclude. 



2 Reduction to Simplified Model 

Spectrum 

As we argued in [11], the only superpartners robustly required by naturalness to lie under 
500 GeV are (in EW gauge basis) the iji,qL = (iL,&z») stops and sbottom, and H u ,Hd 
higgsinos (if their mass arises from a [i term). Along with a SM Higgs boson h, we shall 
consider these the only new particles substantially accessible to the 7 — 8 TeV LHC with 
moderate luminosity. 

Colored superpartners have strong production cross-sections, which suggests that we 
focus our searches on them. As can be seen in Table 1, even these strong cross-sections 
peter out for squark masses above 500 GeV, so that the natural regime for the spectrum 
is also our only hope for direct visibility. We will argue that light higgsinos are typically 
a complication in these searches, either mild or major depending on the spectrum, but 
rarely do they present a spectacular new opportunity. For now we simply neglect them, 
but return at the end of this section to better justify this position. 

The gauge-basis squark states are non-trivially related to the mass-eigenstates after 
EW symmetry breaking, due to two effects, the splitting of II from b^, and the mixing of 
II and tR. The first of these effects is given by the sum of F and D term potentials, 



mr — mr fa . (2.1) 

m ~b L + m h m b L + m h 



The tL — t_R mixing arises from a possible SUSY-breaking ^4-term, Ai c R h u q~L- This results in 
mass eigenstates t\ and I2 (with < m: by convention) related to the gauge-eigenstates 
by an angle 8r. Combining these effects, the sbottom is either the middle or the lightest 
of our squarks. We focus on the former case, with a spectrum ti — b — t\. Note also, 
that even though the mass splitting between b and t\ is essentially a free parameter, it 
cannot be too big if we play by two rules: (a) we keep the lightest squark heavier than the 
top quark, which would otherwise change the phenomenological possibilities, and (b) keep 
the spectrum natural and production cross-sections for b appreciable by not making it too 
heavy. 
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Couplings 



The central novel interaction being considered is the BNV RPV coupling, 



CD 



(2.2) 



enabling stop decay to SM quarks. The exact flavor structure of these couplings is con- 
strained by a variety of low-energy flavor and precision tests. From the viewpoint of LHC 
visibility the central issue is whether the quarks in the dominant BNV couplings carry 
heavy flavor (5r quarks) or not. We have argued in Ref. [11] that in either case, low-energy 
constraints can be satisfied within quite plausible UV flavor paradigms. Relatedly, in the 
case where all three generations of squarks are present, suitable flavor paradigms have also 
been studied by Refs. [30, 31]. We shall therefore consider two cases for stop decay to 
quarks via BNV: (a) one b quark and one light quark, and (b) two light quarks. 

We will always consider the generic possibility that there is at least modest non- 
vanishing A-induced mixing, such that both t± 2 inherit BNV couplings to quarks via their 
tR component. Similarly, they both have weak couplings to the sbottom via their II 
component. Furthermore, mixing also leads to f^t%Z and t^t%h couplings, determined by 
the A-terms and mixing angles. We assume small BNV couplings, X'lu < 1, so that in 
general squarks will only decay through such interactions if decay by W, Z or h emission 
is kinematically suppressed, as for example is obviously the case for t\ (which, recall, we 
are considering as the lightest superpartner). 

Since X'^jj <C 1 is the most straightforward way to comply with low-energy constraints, 1 
it is important to ask how small these couplings can be without resulting in displaced 
vertices at the LHC from a long t\ lifetime. In order to not have a displaced vertex, we 
need c/T to be less than about 1 mm. The expression for the distance traveled before 
decay by a pure tn particle (ignoring mixing) is 



couplings are smaller than this bound, we will have events with jets emerging from displaced 
vertices, which can further help discriminate against background. We will tackle the more 
challenging case in this paper, by assuming that the BNV couplings are strong enough that 
ti decays are prompt. 

Higgsinos 

We are now in a position to understand how higgsinos might affect the LHC physics. We 
have discussed how the spectrum of squarks can be produced and then cascade decay by 

1 Although one can easily satisfy low-energy constraints with X'zij "C 1, there is still a concern that BNV 
can wash out the cosmological baryon asymmetry, if this is generated in the early Universe. One can try 
to turn this into a mechanism for actually generating the baryon asymmetry below the EW scale [32-35]. 
A new robust approach will appear in [36]. 




(2.3) 



Thus, for 300 GeV squarks, we need A" to be roug hly bigger than 2.5 • 10~ 7 . If the BNV 
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EW boson emission, with a final prompt BNV decay to quarks. Higgsinos of comparable 
mass to the squarks allow these steps to potentially be bypassed, by opening up alternative 
squark decays to higgsinos. 

The simplest case would be if the higgsinos were even a little heavier than the stops and 
sbottom. Since direct EW production has substantially lower cross-section, such higgsinos 
would be phenomenologically irrelevant. But if the higgsinos are lighter than the heaviest 
stop, then t 2 decays via EW emission or BNV can be substantially degraded by decay to 
H + b. (The alternate decay to H°t is likely to be phase-space suppressed.) In turn, H + 
will decay (via t\ and BNV) to three jets. In this way, the higgsinos will degrade events 
with leptons from (possibly off-shell) W, Z, and add events with extra b jets. This is the 
basic complication we alluded to earlier: higgsinos can force us to look in multi-jet events, 
without spectacularly high pr, with resonances obscured by combinatorial background, 
and with only the handle of several b jets. 

But fortunately, the higgsinos can easily not degrade sbottom decays even if they 
happen to be lighter than the sbottom but heavier than the lightest stop, because the only 
sbottom decay to higgsinos (for small bottom Yukawa coupling) is to H~t, which is likely 
highly phase-space suppressed. Unfortunately if the Higgsinos are at the bottom of the 
spectrum they will be produced in abundance in t\ — > bH + decays. This decay mode does 
not affect the W production but complicates the resonance reconstruction from the jets. 
This is one of the reasons that we focus on sbottom charged current decays in this paper: 
the phenomenology of sbottom — > W^*H\ is largely "immune" to higgsinos if they are not 
the lightest SUSY particles. 

We proceed by dropping higgsinos from the discussion as part of arriving at our sim- 
plified model of stops and sbottoms in Sections 3 and 4. As discussed above, this will only 
modestly affect our central channel of sbottom production and cascade decay, and is the 
best case for the other channels (but dependent on the spectrum). We will return however 
to the possibility of Higgsinos in the spectrum in Sec. 5 since higgsinos can easily dra- 
matically reduce the contributions of the neutral current decay t 2 — > Z^t±. In particular 
the higgsinos can suppress a yield of neutral current decays in an otherwise spectacular 
multilepton channel. 

3 Signals and strategy of search 

One finds the highest production cross sections for the lightest particles, which would imply 
in our case a search for pair-production of the lightest stop with subsequent decays into four 
jets. However a search for resonances in 4-jet events is very challenging at the LHC [37], 
because it has to deal with a big uncertain QCD background, and even the multijet trigger 
is probably not 100% efficient in this case. 2 Therefore it is fruitful to concentrate on longer 
cascades, which involve Higgs or EW boson emission (either on- or off-shell). This means 
we consider production of heavier states, b, t 2 and subsequent transitions 

b -> W^h, h -> Z(*U U t 2 -> h { *\ . (3.1) 

2 An analogous search for resonances in 6-jet events [38] has a better reach, but it is relevant only for 
light gluinos in an RPV spectrum. 



-6- 



m = 230 GcV 


m = 250 GcV 


m = 270 GcV 


m = 300 GcV 


m = 400 GcV 


m = 500 GcV 


8.3 pb 


5.3 pb 


3.4 pb 


1.9 pb 


0.34 pb 


0.08 pb 



Table 1. Pair production cross sections for stop at NLO for center of mass energy y/s = 8 TeV. 
Sbottom production cross-sections are very similar, slightly bigger though due to electro-weak 
effects. 



While the first process in Eq. (3.1) is fairly robust, the branching ratio of the two other 
processes is model-dependent. The relative rate between the second and the third process 
in Eq. (3.1) is determined by the couplings of the stops to Z and h and by phase space 
effects. While the neutral current decays can have spectacular multi-lepton signature (see 
Sec. 5 for a detailed discussion) it might also happen that the second stop's mass is between 
400 and 500 GeV, rendering the production cross-section tiny (see Table 1). Moreover most 
of the spectacular signatures come from the decays into Z rather than the Higgs, though 
it would be very nice to eventually observe the Higgs in these new physics processes. We 
might however find ourselves in the situation that the mixing angle between the stops, % 
is large and Higgs transitions are preferred. Therefore, it is fair to say that the charged 
current transition is the robust and spectacular channel at the 7-8 TeV LHC and we will 
give it most of our attention. 

Before we continue with a detailed analysis of the cascade decays, we note that 
the table 1 pair-production cross sections for stops were calculated at the NLO with 
Prospino 2.1 [39]. Sbottom production cross sections are usually slightly bigger due 
to electroweak corrections. These numbers will further help us in our numerical estimates. 

Charged-current b — > Wt transitions with a subsequent RPV decay of the stop into 
two jets were first addressed in [25] in the case of resonant production of sbottom. (Pair- 
production was also briefly considered, but was not the primary focus.) Although resonant 
production is not categorically excluded by the bounds on RPV, it requires strong enough 
BNV to raise FCNC and n — n oscillation concerns. However we focus on pair production 
with RPV decay mediated by couplings which can be much smaller than one, and therefore 
safer from low-energy tests. 

In this case the most spectacular signature shows up when both Ws decay leptonically, 
leading to a signature l + l~ jjjj + $t-> where the jets reconstruct two resonances with equal 
masses (see diagram in Fig 1). What should be our search strategy for these events? 
Performing cut-and-count search on events which reconstruct resonances is probably not 
ideal. However we can try to reconstruct resonances with the following steps: 

• Find events with 2 isolated leptons and moderate $t (the latter should be non-zero 
to remove the background from DY dilepton production). 

• Cluster the jets with sufficiently big radius (otherwise there is a danger that we lose 
the hadronic activity which reconstructs the resonance and thereby get edges instead 
of peaks). 

• If the event contains 4 jets (or more), try all possible pairings between the jets, and 
pick up the combination which minimizes the difference between the reconstructed 
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Figure 1. The sbottoms are pair-produced and undergo charged-current decay. When both Ws 
(either on- or off-shell) decay leptonically, they leave a spectacular signature of two leptons + jets, 
which reconstruct two equal-mass resonances. We analyze this signal in Sec. 3 and 4. 

invariant masses. Discard the event if the minimal possible mass difference is too big. 
This step is essentially identical to the standard multi-jet resonances search [40]. 

Unfortunately our events with 2 leptons, MET and multijets have an appreciable back- 
ground, on top of which we are looking for our bumps. This background is heavily domi- 
nated by dileptonic ti (including It\ decay modes). One can show that with an adequate 
choice of cuts all other backgrounds (Z — > t\T\ + jets, DY dileptonic production with jets, 
WW + jets) are highly subdominant to ti, and we will discuss it in more detail in the 
next section. Production cross section for dileptonic ti exceeds our signal by two orders of 
magnitude, and even though the extra jets in these events do not come from resonances, 
reconstructing "by accident" two pairs of jets with similar invariant masses is common. 
The above mentioned steps, plus standard cuts for the overall hardness of the event, are 
still not enough in order to see clear bumps on top of this continuous ti background after 
y/s = 8 TeV run. We therefore use other, less standard discriminators to distinguish the 
signal from the background. 

There are two additional important features which distinguish our signal from the 
background. Usually in a dileptonic ti event, hardness of the entire event correlates with 
the hardness of the leptons and the This happens because the W is often boosted in 
the rest frame of the decaying top. However it is not the case in the signal. As we have 
explained in Sec. 2, naturalness and visibility motivate mild splittings between the stop 
and the sbottom, usually so small that they do not allow emission of the on-shell W. Even 
if emission of the on-shell W is allowed it typically has little boost in the rest frame of the 
decaying sbottom. This results in relatively small pt(0 and $x even if the event overall 
is very hard. We demonstrate the distribution of $t and the transverse momentum of 
the leading lepton in signal and background events on Fig. 2. This immediately suggest 
that just cutting on the tail of high ^j, and high pr{h) should be a decent discriminator 
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P T of the leading lepton in signal and background samples | | E?' 5S distribution in signal and background 




P,(l,) E, mi ' 



Figure 2. Leading lepton pt and $ T distribution in signal and background events. Light red and 
dark red curve stand for the signal spectra 1 and 2 (see Table 2). The blue line represents the 
distribution in dileptonic tt background and the violet line represents the Wi background (which we 
simulate separately since it has slightly different kinematics). See Sec. 4 for details of simulations. 



between the signal and the background. We checked it explicitly and it indeed removes a 
fair portion of the background. We will use a refined version of this discriminator below. 

It turns out one can do even better than just cutting on a high $ T and high pt(Ii) 
events. As we explained, the key feature of the tt events is that usually the leptons and 
the $t are correlated with the hardness of the event, or St defined as 

St = Y,Prtii) + X)pr&) + t'T ■ (3-2) 

i k 

On the other hand in the signal events these quantities are mostly uncorrelated. For this 
purpose we define the following variables: 

_ Pr(h) _ $t /„ Q \ 

n = ^r, r $T = ^- (3.3) 

One should also prefer using these variables rather than pti Pt(Ii) because they are 
dimensionless and therefore cutting on them we do not introduce an explicit scale to the 
problem. We expect these quantities in the signal events to be in general small. We plot 
these variables for signal and background events in Fig. 3 and it follows this expectation. 
Moreover, we see that 77 and r^ T are slightly less dependent on the particular spectrum 
than pt(Ii) and $ T - In the next section we show that using this strategy together with the 
cuts on variables (3.3) we will have an excellent reach after the yfs = 8 TeV, C = 20 fb -1 
run. 

To summarize, the dileptonic channel is an excellent channel for the charged current 
decays. We will elaborate on a feasibility of this search explicitly and make more comments 
on the background behavior and shapes in section 4. 

Finally we briefly comment on semileptonic and all-hadronic decay modes. The latter 
will probably be very hard to utilize, since it just results in multijets (up to 8 or even more) 
events without any evident handles like isolated leptons or missing Et- In the semileptonic 
search one has signal events with isolated lepton, moderate $ T an d at least 6 jets, typically 
resulting in small Pt-i high Ht events. We will not try to elaborate on the feasibility of 
the cut-and-count search in this channel, because evidently these searches are not optimal 
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p T <0<s, " ' E ™";s T ' 

Figure 3. Distribution of n (on the left) and r@ (on the right) as defined in Eq. (3.3). Light red 
and dark red curves stand for the signal spectra 1 and 2 (see Table 2), blue curve represents the 
dilcptonic background and violet represents Iti background. See Sec. 4 for details of simulations. 





m l 


m h 


u{bb*) 


1 


250 GeV 


186 GeV 


5.7 pb 


2 


270 GeV 


189 GeV 


3.7 pb 


3 


300 GeV 


217 GeV 


2.0 pb 



Table 2. Benchmark points for the charged-current decay search. The production cross sections 
are given only for sbottoms . 

(they cannot reconstruct the resonances in jetty channels, taking advantage of the most 
interesting feature of the RPV signal) and basically already exist in some form both in 
ATLAS and CMS collaborations (they do not yield any interesting bounds though). It 
would be interesting to see though how the variables (3.3) can be used in these searches to 
improve further the reach and suppress the backgrounds. 

4 Details of event simulations, backgrounds and reach 

To estimate the feasibility in dileptonic channel as explained in section 3 we analyzed three 
benchmark points points with masses presented in Table 2. For the signal we simulated 
parton-level events with MadGraph 5 [41] for three signal benchmark points given in Table 2 
and showered and hadronized them with Pythia 8 [42]. We wrote down a tailored model 
in Feynrules [43] for MadGraph 5 to capture the effects of the simplified model described 
in Sec. 2. For the background we simulated the events in MadGraph 5 and showered with 
Pythia 6 [44]. In order to capture correctly the effects of extra-jets (which are crucial for 
our analysis) we matched our samples up to two jets with the MLM procedure at 55 GeV. 
Events were clustered with Fast Jet 3 [45, 46]. 3 

Following the discussion in section 3 we reconstruct our events and impose the following 
cuts: 

1. Cluster all the hadronic activity with anti-fey algorithm, clustering radius R = 0.7. 
Relatively large clustering radius is dictated by the fact that we are looking for the 

3 Detector effects are neglected, but the results are sharp enough to survive full treatment. 
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Dijet invariant mass | Dijet invariant mass 




Figure 4. Signal and background events for the benchmark point 1 after C = 20 fb" 1 . Red 
represents the signal, blue the dileptonic tt background, violet is tb, Iti background and grey is 
ti, Tin background. On the LH side plot we do not impose b-veto, while on the RH side plot we 
do. We conservatively assume b-tag efficiency ~ 40%. 

resonances, and smaller radius usually leads to losing relevant hadronic activity. The 
clustering radius is not optimized, but radii of order R ~ 1.0 are likely to be the most 
adequate. 

2. Demand precisely two isolated leptons (carrying more than 85% of the px in the cone 
around the lepton with radius R = 0.3) in each event. We demand pt{Ii) > 20 GeV 
and PT{h) > 10 GeV. 4 The leptons should have \r]\ < 2.5. We discard the event if the 
leptons have same flavor and 81 GeV < mu < 101 GeV to remove the background 
from Z + jets events. 

3. Demand that the event is sufficiently hard, St > 400 GeV as defined in Eq. (3.2) 
and $ T > 35 GeV. 

4. Require four or more hard jets in the event with pt{ja) > 30 GeV. This requirement 
is natural since we are trying to reconstruct two resonances of t±, which both decay 
into two quarks. 

5. Using the variables in Eq. (3.3), demand rg < 0.15 and r\ < 0.15. 

6. Try all possible pairings between four leading jets, and pick up the combination 
which minimizes the difference between the reconstructed invariant masses. Discard 
the event if the minimal possible mass difference is bigger than 10 GeV. 5 If the event 

4 The logic of the cut on the pr of these leptons is dictated by trigger demands. Unfortunately the 
trigger information is not public. However relying on the logic of yfs = 7 TeV run, we hope that the events 
with these leptons should be triggered on with sufficiently high efficiency, namely more than 90% [29]. 
Parenthetically we notice that if the threshold on the px of the leading lepton can be lowered, the results 
that we performed can be further improved. Moreover, some of the events can be triggered on because they 
have sufficient Ht or 4 or more sufficiently high-pr jets. We do not try to take into account the events 
which do not pass these lepton requirement, however lots of them can be "salvaged" since they pass other 
triggers and the ideal search will have to combine several different triggers. 

5 These cuts are not optimized, but it is also not very different from 7.5% of the resonance mass which 
was used in [37] . We explicitly checked our results with respect to variation of this cut. The results are 
rather stable as long as this cut does not exceed ~ 25 — 30 GeV. We leave further optimization of these 
cuts to the experimentalists as it is also going to be affected by jet energy resolution. 
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Dijet invariant mass ~j Dijet invariant mass 




Figure 5. Signal and background events for the benchmark points 2 (up) and 3 (down) after 
C = 20 fb _1 . Red represents the signal, blue the dilcptonic ti background, violet is ti, Iti background 
and grey is ti, t\ti background. On the LH side plot we do not impose b-veto, while on the RH 
side plot we do. We conservatively assume b-tag efficiency ~ 40%. 

has five or more jets with > 25 GeV, try all possible pairings of two and three 
jets. If we get better results when taking the fifth jet into account, use the best 
combination which minimizes the mass difference between the reconstructed objects. 

7. Look for resonances in the reconstructed dijet invariant mass. 

Before we present the results of our simulations we discuss the backgrounds to our 
analysis. Clearly the most formidable background is dileptonic ti (also including the leptons 
coming from leptonic r decays). Naively, one could also worry about Z — > (r/T/) + jets, 
as well as DY l + l~ production and W + W~+ jets. We do not simulate these backgrounds 
and we rely on experimental results which found these backgrounds negligible to ti with 
the cuts which were very similar to ours. First, it was shown in [47] that the background 
Z —> (titi) + jets becomes completely negligible to ti when the third hard jet is required 
in the event. We also see from Table 1 in [47] that the DY background is subdominant to 
ti at least by factor of 5 after requiring at least two hard jets and > 35 GeV. However 
we demand four hard jets in our events, which is supposed to decimate the DY dileptonic 
production and render it completely negligible to ti. Therefore we will further concentrate 
on ti + jets as the dominant background to our signal, and neglect the subdominant 
channels. 

Since we are looking for bumps in the dijets invariant-mass distribution, it would first 
be helpful to understand what effects our cuts have on the backgrounds and how they shape 
the background distribution. Not surprisingly, before all the cuts rrijj in the background 
is a smoothly falling distribution which is peaked around 50 GeV (this peak is carved by 
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our demand from each jet to have pt > 25 GeV. Further demands on hardness of the 
event move this peak to significantly higher values of masses. For example a cut on Ht = 
J2iPT(ji) > 400 GeV (which does a reasonable job with suppressing the backgrounds) 
moves this peak to the vicinity of 200 GeV, which is uncomfortably close to the mass 
scale where we are looking for our resonances. Cuts on the St and the pr of the softest 
necessary jet have a similar effect, therefore we choose the cut in points (3) and (4) discussed 
above to be relatively moderate (one could choose way harder cut which would still remove 
more background than signal events, for the price of moving the peak of the background 
distribution to higher masses). On the other hand, cuts on and r; in point (5) do 
not have this effect, they relatively uniformly discriminate against background from all the 
invariant masses (as this distribution have already been shaped by Ht, St and Pt{j&))- 
Note also, that this terrain of models is still relatively unexplored and the new physics can 
hide at very low masses. It is very natural to expect that the lightest stop has a mass 
of 0(200 GeV) or even lighter. Therefore we emphasize that given a choice of cuts, one 
should always prefer harsher cuts on r^ T and rj, preferring as mild as possible cuts on 
"hardness variables", namely $t> Ht, St,PtU)- This approach is ultimately dictated by 
our attempt to avoid carving spurious bumps on the background distributions. 

Armed with this understanding of the background behavior we turn to the actual 
analysis. For this purposes we assume the NLO ti production cross section at y/s = 8 TeV 
to be 205 pb (this result is taken from MC@NLO [48, 49] with default scale choice). 
Practically, the cross section for ti production is going to be slightly bigger, tallying up 
to 230 pb (from NNLL resummation) with O(10%) uncertainty [50, 51]. However we take 
NLO results to be consistent in our estimates comparing the signal, which we know at the 
NLO, to the background. As we will see our results are strong enough, that increasing the 
background moderately without changing the signal cross section by no mean changes our 
conclusions. We present the results of our simulations in Figs. 4 and 5. In both cases we 
present the results with and without b-veto. As explained in Sec. 2 the lightest stop can 
have decay modes which either include a b-jet or not. Since the dominant background is 
ti, one can achieve much better reach if the signal events contain only light quarks. In this 
case we perform b-veto which further reduces the background, leaving the signal intact. 
The plots on the RH side of Figs. 4 and 5 correspond to this picture. On the other hand, if 
the stop decays to a b-quark + a light flavor, we cannot perform the b-veto (LH side plots), 
but the backgrounds are still under very good control and one has a reasonable discovery 
reach for all three benchmark points. 

Finally, we briefly explain why our resonance search in this channel is much more 
efficient than simple cut-and-count searches. One could naively expect that the discrimi- 
nators that we use to separate the signal from the background are sufficient for an easier 
cut-and-count search. This naive expectation is not true. To illustrate this we plot in Fig. 6 
distributions of St and the invariant mass of the two leading jets. We use all the same 
cuts as in our resonance-searching analysis except the step (6). We see in Fig. 6 that using 
a cut-and-count strategy will be, in the best case, extremely challenging and will require 
detailed understanding of the normalization of the background. Therefore we conclude that 
di-resonance search is the optimal strategy here. 
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Figure 6. Distribution of St (left) and invariant mass of two leading jets invariant mass (right) in 
signal and background events. We use the benchmark point (3) for these plots. Only dileptonic ti 
background is plotted. We use all the same cuts as in resonance search except the point (6). Using 
these distributions as discriminator in a cut-and-count search will be at least very challenging. 

5 Brief Comments on Neutral Current Decays 

Until now we were very detailed in describing the searches for charged-current decays. 
However the charged current decays reveal only part of the full picture. As we explained 
in Sec. 2, in certain spectra the charge-current decay will dominate the collider signatures, 
while in other spectra, the neutral current decays will leave the most spectacular signatures 
(see processes (2) and (3) in Eq. 3.1). We will briefly comment on these processes (see Fig. 7 
for a summary of diagrams) in the current section, but we will be less detailed because as we 
will see the most important potential discovery channel (the multilepton channel) is already 
considered by the CMS. Other channels are not expected to be as strong as mutileptons 
and will mostly favor cut-and-count strategy rather than resonance reconstruction due to 
the high multiplicity of jets. 

In neutral-current decay events we have pair-produced heavy stops decaying into their 
light partners emitting two Z's, two Higgses, or one Z and one Higgs. Z's can be emitted 
either on- or off-shell, while Higgs decays are very unlikely to proceed off-shell due to 
very strong bottom Yukawa suppression. Higgs decays are rarely spectacular, the most 
important Higgs decay mode (assuming that its mass is ~ 125 GeV, as current experimental 
hints suggest) is h — > bb and the third important is h — > gg. If the heavy stop decays to the 
light one emitting two additional jets (b-quarks or gluons) we get a very challenging event 
without any obvious handles, i.e. without MET and/or leptons. The second important 
Higgs decay mode h — > WW* (BR bigger than 20%) and the fourth important (h — > rr) are 
more distinctive, but the resulting going rate of the higgs into two leptons (either through 
W or through r) is smaller than the rate of Z. Therefore we will mostly concentrate on 
neutral current decays with the Z in final state. 

Both Z's decaying leptonically can probably be considered a "golden" channel, even 
though the branching fractions are very small, ~ 0.5%. These rare but very spectacular 
events with up to four isolated leptons and lots of energetic jet activity can be probed by 
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Figure 7. Neutral current decay processes which can be relevant for the most minimal spectrum. 

multilepton searches. The backgrounds for these events from the SM are extremely small 
and therefore even observation of few events can be considered discovery (or, alternatively, 
even the benchmark points which are expected to yield very few events in certain channels 
can be excluded). 

To illustrate this point we explicitly compared a yield of three different benchmark 
points with the results of analysis [29]. This CMS analysis is very special because it 
has a very low pr threshold for the leptons (leptons, which are as soft as pt = 8 GeV 
are considered, higher p^s are required for trigger leptons though). The expected yields 
of all three different points are presented in Table 3. We also compare the yields with 
the theoretical expectations and experimental results, which are quoted from [29]. In all 
these points we assumed, for simplicity of estimation, 100% branching ratio for ti — > t\Z 
decay. However, more realistically, such decays would account for only an order one (but 
highly spectrum-dependent) fraction of ti decays. If higgsinos lie between the two stops, 
than i>2 decay will likely be dominated by t2 —> bH + (mediated by top Yukawa coupling). 
Moreover, as discussed in section 2, the sbottom should lie between the stops, so that ti 
should also undergo charged current decays to sbottom. For squeezed enough spectrum 
such decays may produce even softer leptons than the neutral current decays, and could be 
more challenging to detect. For the estimates below, we ignore such subtleties and study 
the ansatz of neutral-current dominance. We will see that with enough statistics, even with 
significant depletion to other channels, neutral current decays might provide spectacular 
evidence for supersymmetry. 

Clearly with this assumption of 100% branching fraction to Z the first two benchmark 
points in Table 3 are excluded. However full exclusion plots are beyond the scope of 
this paper, because they would demand more refined simulations and more well defined 
assumptions about the Higgsinos, sbottoms and the mixing angle between the stops. Note 
that if the mass splitting between the stops is smaller than 125 GeV the decay will almost 
never proceed through h* due to smallness of bottom Yukawa. Alternatively, pushing the 
heavy stop mass above 300 GeV, we begin to lose the sensitivity in this channel, therefore 
we find that there is a big part of parameter space which is far from exclusion and probably 
is not expected to make a significant contribution to the multilepton channel even after 
£ = 20 ftr 1 run. 

The yields are still smaller than one would naively expect. Most of the events either 
do not have an isolated leading lepton harder than 20 GeV, or the next to leading lepton 
harder than 10 GeV, which renders them unsuitable for a dileptonic trigger. As expected 
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the most populated bins are those with low $t and high Ht-, the Ht in our events comes 
from t% —> jj decays, while the MET is merely instrumental. The prediction in 41, high 
Ht low MET channel are already in tension with the experiment because of extremely low 
(essentially non-existent) SM background. 

It is also expected that there is relatively high chance to lose at least one of the leptons 
due to isolation criteria, or due to high rapidity, so the bins with three leptons turn out 
to be even more informative than the bins with four isolated leptons. Interestingly, the 
only bin where we could predict a significant excess, low MET and high Ht without Z, is 
precisely the bin where CMS observes a non-negligible excess of events, recording 11 events 
where 4.5± 1.5 events are expected. Again, the yield of our benchmark points is too high to 
explain the excess, unless non- vanishing (but also not overwhelming) BR for ti — > H + b is 
assumed. We do not try to claim that it is an anomaly, or that we try to explain this excess, 
however it is probably an interesting channel to watch when i/s = 8 TeV data is analyzed. 
We do not try to estimate the reach of the multilepton channels for i/s = 8 TeV since the 
event yield both of the signal and the background is very low, and the backgrounds are 
very hard to estimate. Nonetheless it is clear that the LHC right now is on the edge of 
probing an interesting region, and high Ht low MET channels are of particular interest to 
watch. 

The second important channel in this category, is one where one of the Z's decays 
leptonically, while its counterpart decays invisibly. This channel has higher branching ratio 
of order 3% but has bigger backgrounds. It also has signatures which naively resemble Im- 
parity conserving SUSY - namely opposite-sign dileptons with jets and Pt- If the mass gap 
between the stops is sufficiently big, such that Z decays on-shell, the signature is leptonic 
Z + jets + pTi naively resembling one of the well-known signatures of R-parity conserving 
gauge-mediation [52]. Alternatively, if the Z decays off-shell we find opposite-sign same- 
flavor pairs and $t, again very similar to a standard R-parity conserving signature with 
decay chain proceeding through a low mass slepton. 

Unfortunately these resemblances are not close enough to be useful. For example, we 
explicitly checked the event yield for all three reference points in table 3 in analysis [53] 
and found that the yields are far below values which one needs in order to have exclusion 
(usually yielding one event or even less in each of the signal regions). These channels are 
so different because they typically have very low missing Et and relatively soft leptons, 
which makes the discovery very difficult with simple cut-and-count experiments. 

Nonetheless one can exploit these events if different strategies are used. If the Z goes 
off-shell we get precisely a signature which is identical to the charge-current decay in the 
dileptonic channel (see Fig. 1) and was analyzed in details in Sees. 3 and 4. One can use 
precisely the same techniques and if an excess in multilepton is found and confirmed in 
the dilepton + MET channel, this can be an excellent cross-check to establish if the excess 
indeed comes from t<i — > Z*t\ decay chains. 

One can also suggest similar searches when Z decays on-shell. In this case the search 
should be modified, because removing the events in the Z-window would wash our signal 
out. Maybe this search is also feasible, however with theory tools it will be hard to 
estimate reliably the background which comes from (Z — > jets with instrumental 
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Table 3. Expected yields of events in the "golden" channel in the multilepton search of [29], 
y/s = 7 TcV with assumption BR(t 2 -> Z*h) = 100%. Channels with high H T , low MET arc the 
most informative. All possible leptonic decays of Z* have been simulated including leptonic rs. We 
define the Z window such that the invariant mass of the OSSF pair is 76 GeV < ran < 106 GeV. 
We do not simulate channels with hadronic rs due to difficulties to mimic one-prong detection 
with our theoretical tools. Three right columns cite the results of [29] where Exp. stands for the 
expected yield of the SM, Err. is the systematic error as it was estimated by the experimentalists, 
and Obs. stands for the observed number of events at C = 4.98 fb . 

Therefore, we point out that this search can be tried, but we do not make any 
conclusions about the backgrounds. 

Although all other channels of Z-decays have much bigger branching ratios, we do not 
see any clear strategies for how these can be utilized. The case where one Z decays lepton- 
ically and the second hadronically would suffer from an enormous Z + Z~ DY background (if 
Z is on-shell, we get Z-\- jets, which is even worse), without even modest MET. One faces 
a similar problem if one of the Z's decays invisibly and the second one decays hadronically. 

6 Conclusions and Outlook 

The main results of this paper are new searches that we propose for the 8 TeV LHC and 
the novel techniques that we find useful to discriminate the new physics signal from the 
background. These searches are motivated by natural SUSY with renormalizable baryon- 
number violating RPV interactions. This yields a set experimental signatures which are 
not efficiently captured by current LHC searches. We pointed out that the signatures of 
charged-current decays can be discovered by a di-resonance search with two additional 
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leptons. The leptons, being soft and accompanied by modest MET are essentially useless 
for cut-and-count search, but provide us an excellent handle and allow us to see the dijet 
resonances despite small production cross section. It is in fact surprising how efficient 
these searches can be. Moreover, there is a good reason to believe that one can do even 
better than our estimates. Although we tried to choose an adequate clustering radius, 
we neither optimized it nor used "grooming" techniques. Simple optimization and using 
"trimming" [54], which is the most adequate grooming technique for these purposes, can 
further improve the sensitivity. 

To efficiently discriminate the signal from the background (which is almost completely 
composed of ti events) we propose to use a set of rather novel cuts combined with more 
standard tools. On one hand we are cutting on the hardness of the MET, leptons and the 
entire event, which is a standard tool, but we also emphasize that these cuts should not 
be too harsh. On the other hand we propose to put an upper cut on ri and rg variables, 
which is a novel way to discriminate signal events where the hardness of the special objects 
in the event (leptons and $t) is uncorrelated with the overall hardness of the event. 

The discriminators that we propose to use are not completely unknown, for example 
CMS was using a much weaker version of 77 as one of the variables in artificial neural 
network in their dileptonic analysis [55] (whose reach is hard to estimate because it uses 
a cumbersome multi-variate approach). The variable r^ T has not yet been use in any 
analysis that we are aware of. We point out that use of these tools can go much beyond 
the particular analyses that we propose, and can be used in cut-and-count experiments as 
well as in resonance searches. We point out that these techniques are suitable in any new 
physics scenario where one finds transitions between the states with small mass splitting 
(see e.g. [56, 57]). 

Finally we point out that the searches that we propose form one more important step 
in the program to map the collider signatures of RPV natural SUSY (for previous works 
see [23, 25]). This is in general a challenging subject, and even R-parity conserving signa- 
tures often demand non-standard approaches [58-62], because regular jets+MET searches 
simply fail. The subject of RPV natural SUSY has received little attention thus far, and 
its collider signatures are still largely unexplored (see however searches by CMS where 
very little or no MET in the signal region is required [63, 64]). It would be interesting 
to study more signatures characterizing natural SUSY with baryon-number violation or 
lepton-number violation, as well as the more challenging spectrum orderings for the sys- 
tem introduced in this paper (lightest superpartner being sbottom or Higgsino). 

We are very hopeful that searches along the lines described above will soon be per- 
formed at the LHC, and will help further our understanding of the grand hypothesis of 
Naturalness. 
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